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Abstract 

We present a detailed report on the decoherence of quantum states 
of continuous variable systems under the action of a quantum optical 
master equation resulting from the interaction with general Gaussian un- 
correlated environments. The rate of decoherence is quantified by relating 
it to the decay rates of various, complementary measures of the quantum 
nature of a state, such as the purity, some nonclassicality indicators in 
phase space and, for two-mode states, entanglement measures and total 
correlations between the modes. Different sets of physically relevant ini- 
tial configurations are considered, including one- and two-mode Gaussian 
states, number states, and coherent superpositions. Our analysis shows 
that, generally, the use of initially squeezed configurations does not help 
to preserve the coherence of Gaussian states, whereas it can be effective 
in protecting coherent superpositions of both number states and Gaussian 
wave packets. 



1 Introduction 

Beyond their fundamental interest in the physics of elementary particles (quan- 
tum electrodynamics and its standard- model generalizations), in quantum op- 
tics, and in condensed matter theory, continuous variable systems arc beginning 
to play an outstanding role in quantum communication and information theory 
Ej , as shown by the first spectacular implementations of deterministic tele- 
portation schemes and quantum key distribution protocols in quantum optical 
settings |S|I1. 

In all such practical instances the information contained in a given quantum 
state of the system, so precious for the realization of any specific task, is con- 
stantly threatened by the unavoidable interaction with the environment. Such 
an interaction entangles the system of interest with the environment, causing 
any amount of information to be scattered and lost in the (infinite) Hilbert 
space of the environment. It is important to remark that this information is 
irreversibly lost, since the degrees of freedom of the environment are out of 
the experimental control. The overall process, corresponding to a non unitary 
evolution of the system, is commonly referred to as decoherence E| ■ It is 
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thus of crucial importance to develop proper methods to quantify the rate of 
decoherence, both for its understanding and for building optimal strategies to 
reduce and/or suppress it. 

In this work we study the decoherence of generic states of continuous variable 
systems whose evolution is ruled by optical master equations in general Gaussian 
uncorrelated environments. The rate of decoherence is quantified by analysing 
the evolution of global entropic measures, of nonclassical indicators and, for 
two-mode states, of entanglement and correlations quantified by the mutual 
information and by the logarithmic negativity. Several initial states of major 
interest are considered. 

The plan of the paper is as follows. In Section [3 we introduce the notation 
and define the systems of interest, together with the quantities wc will adopt 
to quantify decoherence. In Section |21 we introduce and solve the quantum 
optical master equation and its corresponding phase space diffusive equations, 
discussing some general properties of the nonunitary evolution. In Sections I4I7I 
we provide a detailed study of the decoherence of single mode Gaussian states, 
cat-like states, number states and two-mode Gaussian states. Finally, in Section 
|Slwe review and comment the relevant results. 

2 Notation and basic concepts 

The system we address is a canonical infinite dimensional system constituted 
by a set of n 'modes'. Each mode i is described by a pair of canonical conjugate 
operators ii, pi acting on a denumerable Hilbert space Hi- The space Hi is 
spanned by a number basis of eigenstates of the operator fik = aj^a^, 

which represents the Hamiltonian of the non interacting mode. In terms of 
the ladder operators at and one has Xk = (ofe -I- aJ.)/-\/2 and pk = «(aj, — 
ak)/ V^. Let us group together the canonical operators in the vector of operators 
R — {xi,pi, . . . ,Xn,Pn)- The canonical commutation relations regulate the 
commutation properties of the operators: 

[Rk,Ri] = inki , 

where is the symplectic form 

n 
i=l 

The canonical operators Ri may be second quantized bosonic field operators 
or position and momentum operators of a material harmonic oscillator. The 
eigenstates of a; constitutes the important set of coherent states, which is over- 
complete in the Hilbert space Ti^. Coherent states result from applying to the 
vacuum |0) the single-mode Weyl displacement operators Di{a) = e'""*^" 
= A(a)|0). 

The states of the system are the set of positive trace class operators {g] 
on the Hilbert space Ti, = ®'^^{Hi. However, the complete description of any 
quantum state g of such an infinite dimensional system can be provided by one 
of its s-ordered characteristic functions [Jj 

Xs(X)= Tr[gDx]e^"^llV2, (2) 
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with X E R^", II . II standing for the Euchdcan norm R^" and the rt-mode Weyl 
operator defined as 

Dx = e'^^"^, X G ]R2". 

The family of characteristic functions is in turn related, via complex Fourier 
transform, to the quasi-probability distributions Ws , which constitutes another 
set of complete descriptions of the quantum states 

WsiX) ^ \ I d^'^KxsiK) c^^"^^ . (3) 

The vector X belongs to the space F = (R^", fi), which is called phase space in 
analogy with classical Hamiltonian dynamics. As well known, there exist states 
for which the function Ws is not a regular probability distribution for any s, 
because it can in general be singular or assume negative values. Note that the 
value s = — 1 corresponds to the Husimi 'Q-function' W-i{X) = {X\q\X) /tt, 
\X) being a tensor product of coherent states satisfying 

a,\X) = \X) V2 = l,...7i, (4) 

and always yields a regular probability distribution. The case s = correponds 
to the so called Wigncr function, which will be denoted simply by W. Likewise, 
for the sake of simplicity, x will stand for the simmetrically ordered characteristic 
function xo- 

As a meausure of 'nonclassicality' of the quantum state g, the quantity Tg, re- 
ferred to as 'nonclassical depth', has been proposed in Ref. |H] and subsequently 
employed by many authors 




where Sg is the supremum of the set of values {s} for which the quasiprobability 
function Ws associated to the state g can be regarded as a (positive semidefinite 
and non singular) probability distribution. We mention that a nonzero nonclas- 
sical depth has been shown to be a prerequisite for the generation of continuous 
variable entanglement |H| and is strictly related to the efficiency of teleportation 
protocols ^ni- As one should expect, T\n){n\ = 1 for number states (which are 
actually the most deeply quantum and 'less classical' ones), whereas T\a){a\ = 
for coherent states (which are often referred to as 'the most classical' among the 
quantum states) . We note that the nonclassical depth can be interpreted as the 
minimum number of thermal photons which has to be added to a quantum state 
in order to erase all the 'quantum features' of the state. ^ While quite effective, 
the nonclassical depth is not always easily evaluated for relevant quantum states 
(with the major exception of Gaussian states, see the following). 

Therefore, it will be convenient to exploit also another indicator of nonclas- 
sicality, more recently introduced [llj . By virtue of intuition, one should expect 
that remarkable non classical features should show up for quantum states whose 
Wigner functions assume negative values. In fact, for such states, an equiva- 
lent interpretation in terms of classical probabilities and correlations is denied.^ 

^This heuristic statement can be made more rigorous by assuming that a given state owns 
'quantum features' if and only if its P-representation is more singular than a delta function 
(which is the case for coherent states) 

^This is the reason why in the search of CV states able to violate Bell inequalities one is 
lead to consider states with non positive Wigner functions. 
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These considerations have lead to the foUlowing definition of the quantity ^, 
which we will refer to as the 'negative part' of the state g 



^ = j d?^^X\W{X)\-l, (6) 

which simply corresponds the doubled volume of the negative part of the Wigner 
function W associated to g (the normalization of W has been exploited) . 

This work will be partly focused on Gaussian states, defined as the states 
with Gaussian Wigner function or characteristic function x- Such states are 
completely characterized by first and second moments of the quadrature opera- 
tors, respectively embodied by the first moment vector X and by the covariance 
matrix (CM) cr, whose entries are, respectively 



X, = (R,) , (7) 

m{Rj)- (8) 



{RiRj + RjRi) 



The covariance matrix of a physical state has to satisfy the following uncertainty 
relation, reflecting the positivity of the density matrix |12 | 

o- + i^>0. (9) 

The Wigner function of a Gaussian state can be written as 

t^(x) = _^e-5(")"-"'(^-*), ^er, (10) 

ttV Det cr 

corresponding to the following characteristic function 

X{X) = c-5(^--^)^'^(^--*)+»^"^f^^\ . (11) 

A tensor product of coherent states \X) [simultaneous eigenstate of all the 
Oi's according to Eq. is a Gaussian state with covariance matrix cr = ^1 
and first moment vector X. In phase space this amounts to simply displacing 
the Wigner function of the vacuum. 

A single mode of the radiation of frequency lj at thermal equilibrium at tem- 
perature T is described by a Gaussian Wigner function as well. Its covariance 
matrix v is isotrope: v = VI2 with v = [exp(a;/T) -|- l]/[2 exp(w/T) — 2] > 1/2 
(natural units are understood), while its first moments are null. 

The set of operations generated by second order polynomials in the quadra- 
ture operators are especially relevant in dealing with Gaussian states. Such 
operations correspond to symplectic transformations in phase space, i. e. to lin- 
ear transformations preserving the symplectic form VL ^21 ■ Formally, a 2n x 2n 
matrix S correspond to a symplectic transformation (on a 71-mode phase space) 
if and only if 

S'^^S = . 

Simplectic transformations act linearly on first moments and by congruence on 
covariance matrices: cr 1— > crS . Ideal beam splitters and squeezers are de- 
scribed by simplectic transformations. In fact single and two-mode squeezings 
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are described by the operators S^J,r.^ = e^^'^"'"! e'a.aj) ^-^^-^ £ = re'^"^, result- 
ing in single- mode squeezing of mode i for i ~ j ■ Beam splitters arc described 
by the operator Oijfi = e^'^»°^~^'^'°j , corresponding to simplectic rotations in 
phase space. 

A theorem by Williamson |14| ensures that any n-modc CM cr can be written 

as 

cr ^ S^uS , (12) 
where 5 is a (non unique) simplectic transformation and 



■i=i ^ ' 

The Gaussian state with null first moments and CM i/ is a tensor product'^ of 
thermal states with average photon numbers Vi~ 1/2 and density matrices p^^ 



* k=0 \ ' ' 2 
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The set {m;} is referred to as the simplectic spectrum of a, the quantities i/^'s 
being the symplcctic eigenvalues, which are just the eigenvalues of the matrix 
|iri(T|. The uncertainty relation Ineq. (P)) can be simply written in terms of the 
symplcctic eigenvalues 

i^i > 2 = 1, . . . ,n. (15) 

As a last remark about Gaussian states, we briefly address their nonclas- 
sicality. Of course, for the negative part of a Gaussian state one has ^ = 0. 
Remarkably, such an indicator does not detect squeezed states as non classical. 
We point out that this fact is not detrimental to the indicator ^. As a matter 
of facts any Gaussian state can be reproduced in classical stochastic systems 
described by probability distribution, where even an uncertainty relation anal- 
ogous to Ineq. has to be introduced. On the other hand, the nonclassical 
depth of a n-mode Gaussian state g depends only on the smallest (orthogonal, 
not symplcctic) eigenvalue u of the CM cr, which is usually referred to as the 
'generalized squeeze variance' ^3]. The indicator r detects a Gaussian state as 
a nonclassical one (for which r > 0) if a canonical quadrature (possibly resulting 
from the linear combination of the quadratures of the separate modes) exists 
whose variance is below 1/2. The explicit expression for the nonclassical depth 
of a Gaussian state g with CM cr reads 



T„ = max 



1 - 2u 



(16) 



As we have already remarked, coherent states have null nonclassical depth. One 
has to squeeze the covariances to achieve nonclassical features, like subpoisso- 
nian photon number distributions. Regardless of the amount of squeezing, no 
Gaussian state can go beyond the threshold of Tg — 1/2. 

■^as can bo promptly seen from the definition of the characteristic functions, tensor products 
in Hilbert spaces correspond to direct sums in phase spaces 
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In general, the degree of mixedness of a quantum state p of a system with a d- 
dimensional Hilbert space can be characterized by means of the so called purity 
/i = Tt g^, taking the value 1 on pure states (for which = g) and going to l/d 
(that is in infinite dimensional Hilbert spaces) for 'maximally mixed' states. 
The purity is a simple function of the linear entropy Sl = — fJ^d/ {d— 1) and of 
the Renyi '2-entropy' 82 = — In/x, which is endowed with the agreeable feature 
of being additive on tensor product states. While other entropic measures, like 
the Von Neumann entropy, could have been taken into account, the purity has 
the remarkable advantage of being easily computable in terms of the Wigner 
function W{X). Moreover, the global and marginal purities (i.e. the purities of 
the state of the whole system and of the reduced states of the subsystems) have 
been shown to provide essential information about the quantum correlations of 
both two-mode Gaussian states ^1 EI and multipartite, multimode Gaussian 
states ^lElEDl- We also remark that strategies have been proposed to directly 
measure such a quantity, either by quantum networks |23 or by schemes based 
on single photon detections [22]. 

Exploiting the basic properties of the Wigner representation, one has simply 

^i = 'K j w\x) d2"x = ^ / |x(x)p d2"x . (17) 

J 27r J^2^ 

For Gaussian states this integral is straightforwardly evaluated, giving 

^l = — ^ . (18) 

The same result could have been achieved by exploiting Williamson theorem 
and the unitary invariance of This is indeed the way to compute general 
entropic measures of Gaussian states |17| . In particular, the von Neumann 
entropy Sy = —TT:[g\a.g\ of the Gaussian state g is easily expressed in terms 
of the n symplectic Vi 's of the 2n x 2n covariance matrix cr [231 1^ 

n 

Sv = Y.f^^^)^ (19) 

i=l 

with the bosonic entropic function ,f{x) defined by 

f{x) - (x + i) ln(:r +\) - {x \) ln(x - \) . 

This formula will be useful in quantifying the total (quantum plus classical) 
correlations between different modes in two-mode Gaussian states, which will 
be addressed in the following. In general, the total correlations belonging to 
a bipartite quantum state g may be quantified by its mutual information /, 
defined as / = Sv{gi) + Sv{g2) — Sv{g), where gi refers to the reduced state 
obtained by tracing over the variables of the party j 7^ i [^Hl ■ 

Finally, we introduce the definition of logarithmic negativity for bipartite 
quantum states, which will be exploited in the following in quantifying the 
entanglement [i.e. the amount of quantum correlations) of two-mode Gaussian 
states. For such states separability is equivalent to positivity of the partial 
transpose g (PPT criterion) [521 123 The negativity M{g) of the state g is 

*The partial transpose Q is obtained by the bipartite state q by transposing the Hilbert 
space of only one of the two parties. 
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defined as [2Hll2i 

^9) - , (20) 

where ||6|| = Tr V d^d stands for the trace norm of operator 6. The quantity 
M{g), being the modulus of the sum of the negative eigenvalues of g, quantifies 
the extent to which g fails to be positive. The logarithmic negativity E_\f is then 
just defined as E_\f = In || ^|| i. From an operational point of view, the logarithmic 
negativity constitutes an upper bound to the distillable entanglement [2H| and 
is directly related to the entanglement cost under PPT preserving operations 

3 Dissipative evolution in Gaussian environments 

We will consider the dissipative evolution of the infinite dimensional n-mode 
bosonic system coupled to an environment modeled by a continuum of oscil- 
lators. The couplings and the baths interacting with different modes will be 
uncorrelated and generally different, each bath being made up by a different 
continuum of oscillators.. The bath associated to mode i will be labeled by the 
subscript i. The dynamics of the system and of the reservoirs is described by 
the following interaction Hamiltonian 

n „ 

Hint = [wi{Lj)albi{uj)+Wi{Lj)*aibl{uj)]duj, (21) 

4=1 

where bi(uj) stands for the annihilation operator of the ith bath mode labeled 
by the variable uj, whereas Wi{uj) represents the coupling of such a mode to the 
mode i of the system (taking into account the density of environmental modes) . 
The state of the bath is assumed to be stationary. Under the Markovian approx- 
imation, such a coupling gives rise to a time evolution ruled by the following 
master equation (in interaction picture) 

n 

= Y.J + iN^ + 1) L[a,]g - M* D[a,]g + M,D[al]g^ , (22) 

1=1 

where the dot stands for time-derivative, the Lindblad superoperators are de- 
fined as -/^[ojg = 2dgd^ — d^dg— gd^d and D[d]g = 2dgd— ddg— god, the couplings 
are 7, = 27ru'^(0), whereas the coefficients Ni and Mi are defined in terms of 
the correlation functions {bl{0)bi{uj)) — Ni5{u!) and {bi{0)bi{u!)) = Mi5{uj), 
where averages are computed over the state of the bath. The requirement 
of positivity of the density matrix at any given time imposes the constraint 
jMip < Ni{Ni + 1). At thermal equilibrium, i.e. for Mi = 0, Ni coincides with 
the average number of thermal photons in the bath. If Mi 7^ then the bath i 
is said to be 'squeezed', or phase-sensitive, entailing reduced fluctuations in one 
field quadrature. A squeezed reservoir may be modeled as the interaction with 
a bath of oscillators excited in squeezed thermal states several effective 
realization of such reservoirs have been proposed in recent years I34j . In 
particular, in Ref. the authors show that a squeezed environment can be 
obtained, for a mode of the radiation field, by means of feedback schemes relying 
on QND 'intracavity' measurements, capable of affecting the master equation 



7 



of the system More specifically, an effective squeezed reservoir is shown to 
be the result of a continuous homodyne monitoring of a field quadrature, with 
the addition of a feedback driving term, coupling the homodyne output current 
with another field quadrature of the mode. 

In general, the real parameters Ni and the complex parameters Mi allow 
for the description of the most general single-mode Gaussian reservoir, fully 
characterized by its covariance matrix cTiao, given by 



i + iVj + Re Mi Im Mi 



2 



Im M, ^ + N,+ Rc Mi 



(23) 



The non unitary evolution of the single mode system interacting with the reser- 
voir i can be seen as a quantum channel acting on the original state. The 
Gaussian state with null first moments and second moments given by Eq. H23|l 
constitutes the asymptotic state of such a channel irrespective of the initial con- 
dition and, together with the coupling 7;, completely characterizes the channel. 
Now, because of Williamson theorem any centered single mode Gaussian state 
g referring to mode i can be written as 

9 = Sl^.ip.Q'^.Sr,,^, , (24) 

where Sn^ip^ will denote, from now on, the single mode squeezing operator 
Sii^n.Lpi- This fact promptly provides a more suitable parametrization of the 
asymptotic (or 'environmental') state (which is indeed a centered single-mode 
Gaussian state), given by the following equations 

(25) 



^{2N, + 1)2 - 4|M, 



cosh(2r,) = ./1 + 4mL|M,|2, (26) 



tan(2(y5i) 



tan {krgMi 



(27) 



The quantities /iioo, and (fi arc, respectively, the purity, the squeezing pa- 
rameter and the squeezing angle of the squeezed thermal state of the bath. 
The quantity ^ioa is determined, in terms of the parameters of Eq. 1)24(1 . by 
A*ioo = 1/ (2^'i): the purity of a Gaussian state is fully determined by the broad- 
ness of the thermal state providing its normal mode decomposition. 

Eq. (|22ll is equivalent to the following diffusion equation for the characteristic 
function x in terms of the quadrature variables Xi and pi of mode i [7] 



^ 2 

i=l 



[Xi Pi) 



+ [xi pi)cri 



x{x,t). 



(28) 



It is easy to verify that, for any initial condition xq{X), the following expression 
solves Eq. ^ 

x(x,t) = xo(r(i)x)c-^^"'^-(*'^ . (29) 

with the 2n x 2n real matrices F and croo(t) defined as 



m = (g) 
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We mention that Eq. (|22|l can be cquivalcntly recast as a Fokker Planck equation 
for the Wigner function 0. as follows 



Pi J \Op, 



W{X,t). (30) 



Let us now consider a n-mode Gaussian state with CM erg and first moments 
as initial condition in the Gaussian noisy channel. Inserting Eq. in 
Eq. (|29|l shows that the evolving state maintains its Gaussian character and 
is therefore characterized by the action of dissipation on the first and second 
moments. At time t one has 

x{t) ^ r(t)Xo, (31) 
(T{t) = r(t)(Tor(i) + <Too(t). (32) 

In particular, focusing on second moments, Eq. (|32|l is, at any given time t, a 
relevant example of Gaussian completely positive map. Actually, in a more gen- 
eral framework, it can be shown that any evolution resulting from the reduction 
of a symplectic evolution on a larger Hilbert space can be described, in terms 
of second moments, by 

(T ^ X'^aX + Y , (33) 

where X and Y are 2n x 2n real matrices fulfilling Y + iU-iX'^nX > [23EHI- 
Viceversa, any evolution of this kind may be interpreted as the reduction of a 
larger symplectic evolution. 

As a last remark about the dissipative evolution under the master equation 
(|22|l . we point out an interesting general feature concerning a single- mode non- 
squeezed bath, characterized by its asymptotic purity iiao- Let us consider the 
evolution in such a channel of an intial pure non Gaussian state (whose Wigner 
function necessarily takes negative values). It can be shown by a beautiful geo- 
metric argument |39| that the instant tnc at which the state's Wigner function 
gets non negative, so that the nonclassicality of the state quantified by its neg- 
ative part ^ becomes null, does not depend on the chosen state at all. Such a 
time (that is also referred to as 'positive time') reads 

t„e = -fii(l + Moo) . (34) 
7 

In section IHl we will provide a simple proof of this result for an intial number 
state. 



4 Single-mode Gaussian states 

The set of single-mode Gaussian states can be regarded as the simplest contin- 
uous variable arena in which the decay of quantum coherence can be examined. 
The evolution of single-mode Gaussian states in thermal reservoirs has been 
extensively addressed in Ref . [40j , while their dissipative evolution under a gen- 
eral Lindblad equation has been studied in Ref. Ref. [3^1 contains many 
of the results which will be here reviewed for phase-sensitive baths. Both the 
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purity and the nonclassical depth of Gaussian states are completely determined 
by their CM cr, on which we will thus focus. Exploiting again Eq. (|24() . we 
parametrize the 2x2 CM cr through the parameters fi, r and (p, according to 

CTii = ^[cosh(2r) — sinh(2r) cos(2(p)] , 

CT22 = ^[cosh(2r) + sinh(2r)cos(2(p)] , (35) 

(Ti2 = — sinh(2r) su\{2ip) . 

Notice that the purity ^ characterizes the CM according to Eq. (|18|1 . The 
evolution in cL chcLnnel chciractGrizGci by ^, i^qqi ^oo 

and ^pao of an initial state 
parametrized by /io, ''o and (p^) is provided by the single mode (n = 1) instance 
of Eq. (|32|) . Such an equation, together with the paramctrization of Eqs. H35|) 
can be exploited to promptly achieve the time evolution of the parameters r 
and iy9, yielding 



Moo 



(36) 



+ + 2-^ ( cosh(2roo) cosh(2ro) 

Moo ^ 

+ sinh(2roo) sinh(2ro)(cos(2</j^ - 2(^o))) (l - e"-^*) e'^* 

cosh[2r(t)] cosh(2ro) ^ cosh(2r^) , ^. 

7-T = e ^ + (1 - e ^ ) , 37) 

fJ.[t) Mo Moo 

^ _ sinh(2ro) sin(2(^o)e-'^* - sin(2(^^)^ (1 - e"^*) 
tan[2(p(t)] - sinh(2ro) cos(2(^o)e-^* - cos(2^o,)£ (1 - e-/*) ' ^^^^ 

First of all, according to intuition, the purity fi{t) is an increasing function of 
the input purity (Iq: this complies with a general fact about output purities 
of channels, which are maximized by pure states, due to their convexity |38j . 
Moreover, it is immediate to see from Eq. H36() that in a non squeezed thermal 
bath (i.e. for r^o = 0), the purity is maximum at any given time t for ro ~ 0: 
the output purity of such a channel is maximized for rg = 0, that is for a 
coherent input state."' In the theory of measurement, the fact that coherent 
states yield the minimal cntropic production - under non unitary evolution in 
thermal reservoirs ~ is well known and selects such states as privileged 'pointer 
states' in measurement processes |43l I44| .^ 

For phase-sensitive bath, with Too, foo 7^ 0, the purity /i(t) is maximized for 
i~o = ''oo and ipQ = ifoc + 7r/2. This should be expected: in fact, in terms of 
the single mode squeezing operator Sr,ip entering Eq. (|24|l . this means that the 
optimal input state is countersqueezed with respect to the bath, since Sr.ip+f = 
S~^. Indeed, since the purity is invariant under unitary transformation, such a 
result is just a consequence of the fact that the evolution in non squeezed baths 
is optimized by coherent inputs.^ The optimal evolution of purity, plotted in 

^This is a particular instance of a more general result concerning the output purity of 
Gaussian bosonic channels of the form of Eq. 1331 1381 1421 . 

® Notice that the couplings to the bath of oscillators typically considered in these cases are 
not symmetric under the exchange of the two quadratures: this is the reason why, at very 
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Figure 1: Evolution of the purity of various Gaussian states in a channel with fiac = 0.5 and 
rco = 0. The continuous line refers to an initial pure coherent state (^o = 1, = 0), the dashed 
line refers to a squeezed vacuum (^o = 1, ''0 = 1-5) while the dotted line refers to a thermal state 
with ^0 = 0.05 and ro = 0. 



Fig. ^ is simply obtained by inserting tq — Too — in Eq. 

For an initial squeezed input with squeezing parameter tq in a thermal bath 
with Too = (or, more generally, for an initial state with relative squeezing 
^0 ^ foo 7^ 0) the purity fj,(t) may display a local minimum. The condition 
for the appearance of such a minimum can be simply derived by differentiating 
Eq. (|36|l and turns out to be ro > max [/io/^oo, yUoo/Mo]; the time tmm at which 
the minimum is attained can be exactly determined as 
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^-cosh(2ro) 



Mo 



^ - 2cosh(2ro) 



(39) 



The time imm provides a good characterization of the decoherence time of the 
squeezed state: during the initial steep fall of the purity the coherence and the 
information contained in the initial state are irreversibly spread in the environ- 
mental modes. The subsequent revival of the purity is just a result of the driving 
of the state of the system towards the (asymptotically reached) environmental 
one. 

Concerning the nonclassical depth, the smallest eigenvalue u of a single mode 
Gaussian state is simply found in terms of r and ip as u = e~'^^ / (2/x). Insert- 
ing such a result into Eq. Hlt)|) gives the following equation for the non classical 
depth T of a single mode Gaussian state 



T = max 



1 



-,o 



(40) 



small times, some squeezing provides greater purity in such models 1431 . On the other hand, 
the coupling we consider in Eq. I21i is manifcstely symmetric in and pi. 

'^More formally, one can exploit the invariance of the purity under >S'P(2,R) and bring the 
CM (Too of the bath in Williamson standard form: in these canonical basis of phase space the 
channel is non squeezed and coherent states (with CM ctq = I2/2) maximize the purity. To 
go back to the original canonical basis one has to apply the inverse symplectic transformation: 
this explains the previous result about optimization. 
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Figure 2: Evolution of the nonclassicality r in various channels with fj,oo = 0.5. The dotted and 
the continuous lines refer to an initial squeezed vacuum (^o = 1. ''0 = 1) evolving, respectively, 
in a non squeezed channel (dotted line) and in a channel with Too = 0.2 (continuous line). The 
dotted line refers to an initial state with = 0.7 and ro = 1 and the dot-dashed line to an initial 
state with /iq = 1 and ro = 0.5. 



Let us define the quantity K{t) as 

^^^^ ^ cosh(2ro) ^„,, ^ cosh(2r^) _ _ 
Mo Moo 

Notice that n is an increasing function of rg and a decreasing function of /io- 
After some algebra, Eqs. (|37() and (|40|l yield the following result for the exact 
time evolution of the nonclassicality of a single mode Gaussian state 

Tit) = y-^ ^ . (41) 

Such a function increases with both /i(t) and K(t). The choice of the input 
phase of the squeezing which maximizes T{t) at any time is again ipo = ipoo + 
7r/2, maximizing the purity. The maximization of r(<) in terms of the other 
parameters of the initial state is the result of the competition of two different 
effects. Let us consider 7'o: on the one hand a squeezing parameter tq matching 
the squeezing Too maximizes the purity thus delaying the decrease of T{t); on the 
other hand, a bigger value of ro obviously yields a greater initial t(0). However 
the numerical analysis, summarized in Fig. 12 unambiguously shows that, in 
non squeezed baths, the nonclassical depth increases with increasing squeezing 
rp and purity /ip, as one should expect. 



5 Schrodinger cats 

We consider now the following coherent normalized superposition of single mode 
displaced squeezed states 

^ v/2 + 2cos(0) e-2ll^olP ' ^ ^ 
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where |/?o) = S'l-o.o^XolO): B^nd address its evolution under the master equation 
(I22|l . The choice of a null phase in the operator 5^.0 is just a reference choice 
for phase space rotations. 

This state is a relevant instance of cat-like state, i.e. of coherent superpo- 
sition of pure quantum states, whose macroscopic extension has been invoked 
by Schrodinger to illustrate some of the counterintuitive features of quantum 
mechanics More recently, the seminal proposal by Yurke and Stoler |46| . 
besides spurring a great amount of theoretical work aimed at optimizing the 
generation of cat-like states 021, lead to the experimental realization of mcso- 
scopic (ll^oll — 10) superposition of Gaussian states of the radiation field in 
cavity QED The realization of superpositions of Gaussian motional states 
of trapped particles has been demonstrated as well 0^1 , together with the exper- 
imental investigation of their rates of decoherence (SOI- On the theoretical side, 
many efforts have been done to understand and, possibly, suggest methods to 
control the decoherence of such superpositions )51l 1521 1^ 1541 155j . Furthermore, 
we mention that an accurate analysis, under the 'quantum jump' approach, of 
the decoherence of nonclassical quantum optical states (encompassing both cat- 
like and number states) can be found in Ref. where it is also shown how 
nonclassical states may be the result of proper dissipative evolutions. Most of 
the results here reviewed can be found in Ref. |55j . 

Let us define the matrices R = diag ( e''", e"*"") (corresponding to the action 
of 5^0,0 on the 2-dimensional phase space), and ctq = l/2i?^. The Wigner 
function associated to the state \Po) reads 



1 



471(1 + 003(6*) c-ll-^iill')VDet ctq 

-i(X^-A'f7R)<T-i(X-RXo) , -i(X^-|-Xo^R)<T-i(X+RXo) 



consisting in the two Gaussian peaks at the phase space points Xq and —Xq, 
linked in phase space by the oscillating interference terms. Obviously, this 
Wigner function is non positive. However, formally, such a function is just the 
sum of four displaced Gaussian terms. The linearity of the considered dissipative 
evolution permits to simply solve the evolution of the cat state, by following the 
evolution of its four Gaussian terms according to Eqs. (j^l^ . One gets 



47r(l cos(6l) e-ll^t,ll')v/Detcr(t) 

+ (,-^(x^+o"^'J>^:JR)<T(t)"'(x+c-^'RXo) 

_^ e""-^""' |^g-i(X^-io-^*Xo^a.R)<T(t)-i(X-l-jo-i'Ru;Xo)+ie _|_ ^^^^^ 

(44) 

where (T{t) is given by Eq. H31(l with ctq defined above. 

Figure 13 provides a relevant example of dissipation of a cat state in a ther- 
mal environment, isotrope in phase space. The negative part ^ of the Wigner 
function reaches the value at a time t„c — 0.47"^, in agreement with Eq. (|34|l . 
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Figure 3: Evolution in phase space of the Wigner function of an initial non squeezed cat-like 
state with Xq = (1 1) and 9 = in a thermal channel with /ioo = 0.5 at times t = (a), 
t = 7^^/2 (b), t = (c) and t = 4'y~^ (d). Darker colors stand for lower values, the scale of 
each plot is normalized. The negative lobes (in which the Wigner function takes negative values), 
evident in subfigure (a) are already disappeared in subfigure (b). Actually, the positive time tnc 
of such a reservoir is tnc — 0.47^^ (see Eq. 1341 '). 

As already mentioned, this time is feature of the bath and does not depend on 
the initial pure (non Gaussian) state. 

The exact analytical expression of the purity of the evolving superposition 
is easily determined by Gaussian integrations, according to Eq. p7(l 
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Figure 4: Evolution of the purity of initial cat-like states. The asymptotic purity of the channel 
is = 0.5. The dotted line refers to a cat state with Xq = (1, 1) in a non squeezed channel. 
The dashed and the continuous lines refer to an initial non squeezed cat with = (100, 100) 
evolving in a non squeezed channel and in a channel with Too — 0.88 and (/3oo = —tt/8. The dot- 
dashed line refers to an initial state with A'^ = (10, 10) and ro = 2 evolving in a non squeezed 
channel. 



M/3o.eW = 8(l + cos((?)e-ll^"ll")VDeto-(t) 



2(1 + 6-" '"^oSmXo^ _^2c-2|l'^'of (^008(26*)+ I 



-T"xJ'T(t)Xo 



-h4e" 



11^0 



m (e 



'xJ'jXoTr[JS(t)]*/4 



4- c.c. 



) 



with 



S(0 = Rcr(t)-iR , T(t) = (Dot cr)-iS(t)- 



J = 



1 i 

i -1 



(45) 



(46) 



Eq. (|45|) shows that the decoherence rate increases with the 'dimension' of the 
cat, quantified by ||^oll; in the hmiting instance Xq = 0, Eq. H45() reduces 
to Eq. (|36|l for an initial squeezed vacuum, which decoheres more slowly than 
the equally squeezed cat-like states. Moreover, in general, the terms depend- 
ing on the coherent phase 9 are suppressed by exponential terms of the form 
exp(— ||Xo|p), so that the decoherence rate in terms of the purity is only slightly 
influenced by the choice of 9. Examples of decoherence of cat states can be seen 
in Fig. 21 In all the instances the purity displays a fast initial fall, during which 
all the coherence and the information of the pure cat-like state are lost. The 
typical time scale in which the minimum of the purity is attained is in good 
agreement with the estimate tdec = /'^W^oW'^ ^ holding for the decoherence 
time of a cat state in a thermal bath [SUI ■ As can be shown analytically |55] , the 
phase space direction of the cat, determined by the angle ^0 = arctan(a;o/po)j 
providing the maximal delay of decoherence at short times {i.e. for "ft ~ 1) is 
given by S,o = + ""/S for a squeezed cat in a non squeezed bath or, equiv- 
alently, by ^0 ~ foe for a non squeezed cat in a squeezed bath. These two 
instances are, as already noted, unitarily equivalent. In general, the evolution 
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Figure 5: Comparison between the evolution at short times {i.e. for 7tdec~i) °f the purity of 
an initial non squeezed cat (continuous line) and that of squeezed cats with optimal choice of the 
optical phase. In all instances = (4, 4), ficc = 0.5 and = 0. The dashed line refers to a cat 
state with rg = 1, whereas the dotted line refers to a state with ro = 1.5. The decoherence time 
of such cats can be estimated as t^^,, ~ 0.037^^, in good agreement with the decrease of purtiy 
of the non squeezed cat. The remarkable delay of decoherence induced by squeezing the cat can 
be appreciated, especially at t ~ tdec- 

of the purity of an initial state in a squeezed reservoir is identical to the one of 
the counter-squeezed initial state in a thermal reservoir. Therefore, the same 
protection against decoherence granted by squeezing the bath can be achieved 
by orthogonally squeezing the initial state. Indeed, with the optimal, previ- 
ously discussed, locking of the optical phase, an optimal value of the squeezing 
ro maximizing the purity in non squeezed baths does exist. As illustrated by 
Fig. |5l squeezing the initial cat (or the bath) can provide a significant delay 
of the complete decoherence of the cat state, better preserving the interference 
fringes in phase space. 

6 Number states 

As a last example of single mode state we quantify the decoherence of number 
states |n)(n|. Such states can be considered as probe of fundamental quantum 
mechanical features and are also required in several quantum communications 
tasks I^EHj. Different methods for the generation of Fock states have been pro- 
posed, both for traveling-wave and cavity fields. For traveling-wave fields, these 
methods are principally based on tailored nonlinear interactions |59| , conditional 
measurements [^01, state filtering or state engineering [Hlj. A further pos- 
sibility to generate number states with high fidelities by atom-field interactions 
in high-Q cavities has been recently suggested [HSl- The actual experimental 
generation in quantum optical settings seems to be at hand, by both determin- 
istic I^ESI and probabilistic ('post-selective') schemes (and the techniques 
to realize such states for motional degrees of freedom are well mastered [HZI), 
even if the numerical analysis suggests that environmental decoherence could 
still hamper the very possibility of generating pure number states [BHl- These 
reasons motivated an accurate investigation of the decoherence rate of number 
states, carried out in Ref. [231 ■ We review such results, adding the analysis of 
the nonclassicality of the evolving states. 
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The characteristic function x„ associated to the state is promptly 

found and reads [7| 

XniX) = {n\D^\n) = c-^L„(||Xf ) , (47) 

where L„ is the Laguerre polynomial of order n: Ln{x) ~ 'J2m=o ^^m\ im) ■ 
that, exploiting Eq. I|29(l . one at once finds the evolution of such an initial state 
in the channel 

Xn(0-in(^e--*) e-i^^"-«^, (48) 

with 

'^W = ^e-'^* + '^oo(l- e-^*). (49) 

According to Eq. p7|) one can then determine the purity /i,i(t) of the evolving 
number state [701 

^^,,{t)=e^' r e-^^L„{s)Io(^-^^^^^{e'''^l)s) ds , (50) 

Jo \ ^Moo / 

where Io{x) — Jo{ix) — X^tLo (2^^1)2 the zero order modified Bessel function 
of the first kind and 

^ _ e'y* + /ipo - 1 

For a thermal channel, with Too = 0, such an expression can be further simplified 
to achieve an exact analytical expression for the purity, yielding |7()j 

Mi)=e^*^^P„(l + ^^) , (51) 

where P„ is the Lcgendre polynomial of order n: Pnix) = 2^.Tx^(-'^'^ ~ -'^)"- 
Again, we point out that the squeezing of the bath has the same effect on 
the purity as the counter squeezing of the initial number state, amounting to 
consider a 'squeezed number state'. The numerical analysis of Eq. H5U|I at short 
times (for < I) shows that ^n(i) is a decreasing function of rao- the squeezing 
of the bath does not help to preserve the coherence of number states. Also, the 
purity at any given time is a decreasing function of n: number states of higher 
order are more fragile and decoheres faster. 

Let us now deal with the evolution of the negative part ^ of a number state 
\n), quantifying the dccohcrcncc effect on the nonclassical features of the state. 
The initial value of such a quantity increases with increasing n (higher order 
number states are regarded as 'less classical' by this indicator). Subsequently, 
during the dissipation in the bath, the negative part ^ decreases up to a time tnc 
- determined by Eq. (|34|l - at which it reaches the values and the nonclassical 
features of the state related to ^ are erased. Interestingly, a direct determination 
of the time tnc can be easily provided for the relevant instance of number states 
evolving in non squeezed thermal baths (with Too =0). In such a case, the 
spherically symmetric characteristic function of Eq. H48|l for Too = can be 
Fourier transformed to get the Wigner function Wn (t) 

W^n t = :) I ^1 c~ 



-26-^^11X11 



(52) 
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with 



C(t) = — [l-(l-Moo)e-^*] and = — [l - (1 + Moo) 6"^*] . 

Moo Moo 

Since Laguerre polynomials of any order have positive roots and are always pos- 
itive for negative arguments, Eq. H52(l implies that the time i„c is determined by 
the condition ?/(tnc) ~ 0, yielding i„c = 7^^ lii(l-l-Moo)- This result is just a spe- 
cific instance of Eq. H34() , which can be applied at any pure non Gaussian initial 
state. It can also be found in Ref. [71], where the remarkable independence of 
the time tnc on the order n of the number state had already been stressed. The 
evolution in phase space of the Wigner function of Eq. (|52|l is showed in Fig. 
Fig.Elshows the time dependence of the negative part ^, numerically integrated 
for the first four number states in a thermal reservoir. Even though the initial 
negative part increases with increasing n, the quantity £^{t) is not increasing 
with n at any time: indeed, lower order states better preserve such nonclassical 
features when approaching the time tnc (which, we recall once again, does not 
depend on the initial pure non Gaussian state). 

A relevant instance to exemplify the decoherence of number states is provided 
by the coherent normalized superposition \ipoi) = {\0)+ e*''|l))/-\/2, constituting 
a microscopic Schrodinger cat. The characteristic function xai of this state is 
simply found [J] 



Xoi(a) = [2 - e-^*|a|2 - e-¥(a* e"*'' - ae*")] . (53) 

Inserting xoi ^-s the initial condition in Eq. H29|l and performing the integration 
of Eq. (|17|l yields, for the purity of the initial cat-like state evolving in the 
channel 



Moi 



(i, r) = Aiy- e-^'<'^(iioo + ( e^* - l)(cosh(2r) 
2moo ^ 

+ cos(2??- 2.^)sinh(2r)) 

5 



2-(4ML + 8(e^*-l)MooCOsh(2r) 
2mSo ^ 

+ (e'^*- 1)2(3 cosh(4r) + l)) (54) 



where 



(1 - c-'^*)^ + e"27* + 2— cosh(2r) 

Moo Moo 



1-1/2 

(55) 



is the purity of an initial vacuum in the channel, found in Sec.01 Eq. 1)54(1 shows 
that the evolution of the coherent superposition is sensitive to the phase Lp of 
the bath. It is straightforward to see that the optimal choice maximizing purity 
at any given time is provided hy "d = (p -\- 7r/2. Fixing such a choice, we have 
numerically analyzed the dependence of moi on the squeezing parameter rao- For 
small Too the purity moi increases with r. The optimal choice for r^o depends on 
time, for = 0.5 it turns out to be r ~ 0.28. The relative increase in purity for 
several choices of the squeezing parameter Too is plotted in Fig. |H1 as a function 
of time. It is interesting to compare this analysis of decoherence with the one 
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Figure 6: Evolution in phase space of the Wigner function of the initial number state |2) in a 
thermal channel with /ioo = 0.5 at times t = (a), t = 7^^/4 (b), t = -y^^ (c) and t = 1.57^'^ 
(d). Darker colors stand for lower values, the scale of each plot is normalized. The time t„c, at 
w/hich the Wigner function of this state gets positive is tnc — OA'y~^. As can be seen, att = 7~^, 
the central minimum deriving from the initial negative zone is still evident, but takes only positive 
values. 



previously carried out for Gaussian catlike states. Indeed, notwithstanding the 
deeply quantum nature of a superposition of number states, its decoherence 
rate is comparatively slow. Actually, the purity of the considered superposition 
in a thermal channel reaches the asymptotic value of the channel, after the 
initial decrease, in a time t ~ 0.67"^. Such a time length corresponds to the 
decoherence time tdec of a superposition of two Gaussian terms displaced in 
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Fi gUre 7: Time evolution of the negative part ^ of the number states |1) {dotted line), |2) (dot- 
dashed line), |3) (dashed line) and |4) (continuous line), in a thermal reservoir with /ioo = 0.5. 
For such a reservoir, the Wigner function gets positive at tnc — 0.47^^^. 



phase space of only one coherent photon (in opposite directions with respect to 
the origin, i.e. with ||Xo|P = 1 in the notation of the previous section). Despite 
the relevant intrinsic differences between these two kinds of Schrodinger cat 
states, their decoherence is basically driven by the same process, due to the 
entanglement of the system with the environmental degrees of freedom. 

We remark that the time of decoherence can be much shorter of the time 
characterizing the energy relaxation |31II5()| . which constitutes however a strict 
upper bound on the former. This fact is a manifestation of a general feature of 
quantum mechanics. Nonclassical superpositions decohere on a time scale of the 
photon lifetime in the channel, regardless of the other parameters: once a single 
photon is added or lost, all the information contained in the original state leaks 
out to the environment. This can be understood, euristically, by considering the 
action of the annihilation operators a which, in general, modifies the coherent 
phase of the superposition. Therefore, as soon as the probability of losing a 
photon reaches 0.5, the original superposition turns into an incoherent mixtures 
of states with different phase, whose interference terms cancel out each other 
[561 172| . No coherent behaviour can survive such a dissipative process and be 
afterwards revealed by interferometry. 

7 Two-mode Gaussian states 

Two-mode Gaussian states are the simplest example of continuous variable bi- 
partite states. Their decoherence under the quantum optical master equation 
can be therefore characterized also by investigating the evolution of the cor- 
relations between the two modes of the systems. In particular, the decay of 
quantum correlations, i.e. of the entanglement, quantified by the logarithmic 
negativity, may be adopted as an indicator of decoherence. Due to their clear 
interest, concerning both applications in quantum information and the study 
of fundamental features of entanglement, the behaviour of two- mode Gaussian 
states under non unitary evolutions has attracted a remarkable theoretical in- 
terest in later years I73[731|7Sl|7ni[771[73|7ni|HDl. We review here the results of 
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Figure 8: The relative increase in purity, defined by Afj,/^ = {fj,oi{t,r) — /^oi (t, 0))//ioi (tj 0), 
as a function of time during the evolution of the superposition \4'oi) in Gaussian channels. The 
optimal condition "d = ip + tt/2 is always assumed, while fioo = 0.25. The solid line refers to a 
bath with r = 0.28, close to the optimal value; the dotted line refers to a bath with r = 0.4 and 
the dot— dashed line refers to a bath with r = 0.1. 



Ref. [HOI; moreover, we consider the instance of different couplings to the bath 
and provide a detailed study of the evolving nonclassical depth. 

Before addressing the analysis of their dccohcrcncc in detail, let us recall 
some basic facts about two-mode Gaussian states. The 4x4 covariance matrix 
(T shall be conveniently written in terms of the three 2x2 submatrices a, (3, 7 



a 7 

7^ (3 



(56) 



The CM cr can be put into the so called standard form a-gf through a local 
symplectic operation Si — Si (B S2 



STtS, 



(a Cl \ 

a C2 

Cl 6 

V C2 b j 



(57) 



In what follows, let us suppose |c2| > |ci|. States whose standard form fulfills 
a = b are said to be symmetric. Let us recall that any pure state is symmet- 
ric and fulfills ci = — C2 = ^ 0? — 1/4. The correlations a, 6, ci, and C2 are 
determined by the four local symplectic invariants Deter = (ab — Ci){ab — c|), 
Dcta = a^, Dct/3 = 6^, Dct7 = C1C2. Therefore, the standard form corre- 
sponding to any covariance matrix is unique (up to a common sign flip in the 

Ci's). 

The 5'p(4 R) invariants Deter and A(cr) = Dcta + Dct/3 + 2 Dct7 permit to 
explicitly express Ineq. © in terms of second moments 



A(cr) < - +4 Deter 

4 

and determine the symplectic spectrum {vz^} of er, according to |24| 



(58) 



A(er) T VA(cr)2 -4 Deter 
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A relevant subclass of Gaussian states we will make use of is constituted 
by the two-mode squeezed thermal states. Let Sr = Si2,r,o be the two-mode 
squeezing operator between the modes 1 and 2 with real squeezing parameter r 
and let i/^ be the tensor product of identical thermal states of global purity /x, 
with CM i/p = l/(2y//i)l- Then, for a two-mode squeezed thermal state ^^^^ 
we can write ^^i^^ ~ Sri^^iSl- The CM of ^^^^ is a symmetric standard form 
satisfying 

cosh2r sinh2r 
a = , ci = -C2 = . (59) 

201 

In the instance /i = 1 one recovers the pure two-mode squeezed vacuum states. 
Two-mode squeezed states are endowed with remarkable properties related to 
entanglement [811 . in particular they are the maximally entangled states for 
given marginal and global purities |16l I17| . 

We recall that the necessary and sufficient separability criterion for two-mode 
Gaussian states is positivity of the partially transposed density matrix ("PPT 
criterion") [SHI. It can be easily seen from the definition of W{X) that the 
action of partial transposition amounts, in phase space, to a mirror reflection of 
one of the four canonical variables. In terms of the Sp^^s. invariants, this results 
in changing the invariant A(cr) into A((t) = A(cr) = Det a + Det /3 — 2 Det 7. 
Now, the symplectic eigenvalues Dzp of the partially transposed CM cr read 



\ 



A(<t)t \/A(cr)2 -4Dctcr 

■ (60) 



The PPT criterion then reduces to a simple inequality that must be satisfied by 
the smallest symplectic eigenvalue of the partially transposed state 

^->l^ (61) 

which is equivalent to 

A(cr) < 4Detcr + i . (62) 

The above inequalities imply Det 7 = C1C2 < as a necessary condition for a 
two-mode Gaussian state to be entangled. The quantity i>_ encodes all the 
qualitative characterization of the entanglement for arbitrary (pure or mixed) 
two-modes Gaussian states. Note that i>_ takes a particularly simple form for 
entangled symmetric states, whose standard form has a ~ b 



= v/(a-|ci|)(a-|c2|) . (63) 

The logarithmic negativity Ej^ of two-mode Gaussian states is a simple function 
of D-, which is thus itself an (increasing) entanglement monotone; one has in 
fact [HI 

£■^(0-) = max{0, -ln2£>_} . (64) 

This is a decreasing function of the smallest partially transposed symplectic 
eigenvalue quantifying the amount by which Inequality (jBlj) is violated. 
Thus, for our aims, the eigenvalue i>- completely qualifies and quantifies the 
quantum entanglement of a two-mode Gaussian state cr. 
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The smallest eigenvalue u of crgf (which determines the nonclassical depth 
T according to Eq. (|16|l ) is easily determined 



2u a + 6- -^(a- &)2 +4ci , (65) 

reducing to m = a — |c2| for symmetric states and to u = e^^''/(2y//i) for two- 
mode squeezed thermal states. 

The evolution of two-mode Gaussian states in the noisy channel is described 
by Eq. I)32|l with n = 2. The channel is completely determined by the quantities 
A*ic»i fioo, fioo and 7i, for i = 1,2. Notice that, if 71 7^ 72, then a change in the 
values of the couplings to the bath 7i's does not reduce to a rcscaling of time and 
may significantly affect the evolution of the relevant quantities in the channel. 
For the study of the cntropic measures and of correlations, we will restrict to 
initial states in the standard form of Eq. 1)57(1 , with no loss of generality since all 
such quantities are invariant under local unitary operations. On the other hand, 
the nonclassical depth r is not invariant under such operations. Determining the 
evolution of such a quantity in the general instance is slightly more involved. 
For the sake of simplicity, we will study such evolution in relevant instances, 
which can be conveniently handled and illustrate the general behaviour of the 
nonclassical indicator. Henceforth, we will set ipioo = as a reference choice for 
phase space rotations. 

Exploiting the results we have just reviewed, together with the general def- 
initions of Sec. El we can determine the exact evolution in the channel of the 
entropic measures fi and Sy-, and of the quantum and total correlations, re- 
spectively quantified by E_\f and /. In Appendix we provide the explicit 
expression of the time dependent terms, allowing to compute such evolutions, 
in the instance of equal couplings: 71 = 72 = 7. 

As for the evolution of the purity ^ and of the von Neumann entropy Sy — 
whose decrease quantifies the information which the composite two mode state 
'as a whole' loses by interacting with the environment - some analytical state- 
ments can be done. It can be shown by means of a variational approach |38j that 
the purity of a given channel of the form of Eq. H32|l is maximized by an uncor- 
related state (with Ci = C2 = in our notation). Its maximization is therefore 
achieved by the (obviously separable) product of two 'countcrsqueezcd' states, 
which, as we have seen in Sec. ^maximizes the local purity relative to the two 
single-mode channels.^ The optimal purity evolution reduces therefore to the 
square of the optimal purity evolution for single mode channels, previously stud- 
ied. This feature holds for any value of 71 and 72. An analogous argument can 
be applied to the von Neumann entropy Sy which, we recall, is fully determined 
by the quantity limp_,o Tr gP. However, so far, the fact that the minimal Sy at 
any given time is achieved by an uncorrelated input has been proved only for 
7i =72- The numerical analysis, summarized in Fig. |3 remarkably supports 
the conjecture of the additivity of the minimal output von Neumann entropy 
also for 71 ^72.^ 

We now move to consider the decay of the entanglement between the two 
modes of the field, i.e. the leaking to the environment of the information con- 

^This is a particular instance in which, rstricting to the Gaussian setting, the maximal 
output purity of a tensor product of channels is 'multiplicative' 1381 . 

^The additivity of the minimal von Neumann entropy corresponds to the multiplicativity 
of the maximum of the quantity limp^o Tr gP. 
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Figure 9: Time evolution of the von Neumann entropy in a thermal channel with 71 = 1, 
72 = 2 and ^loo = M2oo = 0.25. The continuous line refers to the conjectured optimal evolution, 
achieved by a pure separable input state with a = b = 1/2 and ci = C2 = in a non squeezed 
channel; the dotted line refers to a n inital pure two mode squeezed state with r = 0.5 in the 
same channel. The dashed and dot-dashed lines refer to a squeezed channel with riao = ''200 = 1 
and, respectively to an initial pure two mode squeezed state with r = 1 and an intial thermal two 
mode squeezed state with = 1/16 (equal to the asymptotic purity). 



taincd in quantum correlations between the two modes. Supposing that the 
couplings to the two baths are equal (71 = 72 = 7) and making use of the 
separability criterion given by Ineq. (|62|l . one finds that an initially entangled 
state becomes separable at a certain time t if 

u e"*"^* + V e-3-^* + w e-^"^* + 2/ e"''* + z = . (66) 

The coefficients u, v, w, y and z are functions of the nine parameters character- 
izing the initial state and the channel (see Add. Eq. (|66|l is an algebraic 
equation of fourth degree in the unknown k = e^''*. The solution fcent of such 
an equation closest to one, and satisfying kgnt < 1 can be found for any given 
initial entangled state. Its knowledge promptly leads to the determination of 
the "entanglement time" tent of the initial state in the channel, defined as the 
time interval after which the initial entangled state becomes separable 

tent = - - In kent ■ (67) 

7 

The entanglement time tent can be easily estimated for symmetric states 
(for which a = 6) evolving in equal thermal baths {i.e. with 71 = 72 = 7 and 
fJ-ioo = A*2oo = y^fJ-oc)- In such a case the initally entangled state maintains its 
symmetric standard form during the time evolution. Recalling that |ci| < |c2|, 
we have that Eqs. H61() and H63() provide the following bounds for the entangle- 
ment time 

<7tent<\n(l + ^?M:^^) . (68) 




^''clearly, in the general instance of different couplings (71 ^ 72), Eq. 1661 would turn in a 
system of fourth degree in the two unlnown e 'i' and e '2'. Such a situation does not pose 
any conceptual problem and can be treated in much the same way as the one here described, 
by explicitly determining the coefficients of the system. 
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Fi gure 10: Time evolution of the logarithmic negativity of an initial two mode squeezed thermal 
state with fi = 0.8 and r = 1 in several channels with 71 = 72 = 7- The continuous line 
refers to a non squeezed bath with /iioo = M2oo = 0.5 (corresponding to 0.5 thermal photons); 
the dotted line correspond to a thermal bath with fiiao = 0.25 and fi2cxi = 1 (with a global 
asymptotic purity equal to the previous one); the dashed and dot-dashed lines refer to a bath with 
Mioo = A'2tx) = 0.5, rioo = r2aa = 1 and 95200 = ((/P2oo = 7i"/4) for the dashed (dot-dashed) 
line. 



Note that is the global purity of the asymptotic two mode state. Imposing 
the additional property ci = — C2 amounts to consider standard forms which can 
be written as squeezed thermal states (see Eas. I59|l. For such states, Inequality 
reduces to 



In particular, for /i = 1, one recovers the entanglement time of a two-mode 
squeezed vacuum state in a thermal channel [271 1771 17^ . We point out that 
two-mode squeezed vacuum states encompass all the possible standard forms 
of pure Gaussian states. 

The results of the numerical analysis of the evolution of the logarithmic 
negativity for several initial states are reported in Figs. 1101 andllll In general, 
one can see that a less mixed environment better preserves entanglement by 
prolonging the entanglement time. More remarkably. Fig. ^| shows that a local 
squeezing of the two uncorrelated channels does not help to preserve the quan- 
tum correlations between the evolving modes. Moreover, as can be seen from 
Fig. ^2 states with greater uncertainties on, say, mode 1 {a > b) better preserves 
its entanglement if bath 1 is more mixed than bath 2 (/iioo < M2oo)- Fig. 1111 
also shows that, even for initial non symmetric states, unbalancing the cou- 
plings to the two single mode reservoirs (while leaving their average unchanged: 
7 = (71 + 72)72) only slightly affects the evolution of the entanglement in the 
channel; an accurate numerical analysis shows that a greater coupling to the 
more mixed initial mode {e.g., 71 > 72 if a > b) enhances the preservation of 
the initial quantum correlations. Also, for symmetric states evolving in squeezed 
baths, one can see that the entanglement of the initial state is better preserved 
if the squeezing of the two channels is balanced. 

An interesting feature concerns the evolution of the mutual information /, 




(69) 
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Figure 11: Time evolution of the logarittnmic negativity of an inital entangled non symmetric 
state, obtained from the squeezed thermal one considered in Fig. llOl bv adding 0.2 to the element 
a of the standard form (added noise on mode 1 quadratures.) The solid line refers to a bath with 
71 = 72 = /^loo = A'2cx) = 1; the dotted line refers to a channel w/ith 71 = 0.5, 72 = 1.5 and 
Mioo = M2tx3 = li the dashed (dot-dashed) line refers to a bath with 71 = 72 = 1, /^loo = 1/9 
(/iloo = 1) and fi2oo = 1 (M2oo = 1/9). The label 7 is defined by 7 = (71 -|-72)/2. 
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Figure 12: Time evolution of the mutual information of Gaussian states in an environment with 
71 = 72 = 7 and fJ,i = fJ.2 = 1/3. The continuous line refers to an entangled state with a = 2, 
b = ci = — C2 = 1 in a non squeezed environment; the dotted line refers to the same state in an 
environment with ri = r2 = 1; the dashed line refers to a non entangled state with a = b = 2, 
ci = — C2 = 1.5 in a non squeezed environment; the dot-dashed line refers to the same state in a 
squeezed environment with ri = r2 = 1. The squeezing angle ip2 has always been set to 0. 



illustrated in Fig. El for some relevant cases: at long times, such a quantity 
is better preserved in squeezed channels. This property has been thoroughly 
tested both on non entangled states, featuring only classical correlations, and 
on highly entangled states, and seems to hold generally. 

The instance of a standard form state in a tensor product of two thermal 
channels (parametrized by 7^ and ^^oo, for i = 1, 2) is especially relevant, since it 
gives a basic description of dissipation in most experimental settings, like fiber- 
mediated communication protocols. A simple analysis straightforwardly shows 
that in this instance both the purity and the logarithmic negativity (that is, the 
entanglement) of the evolving state are increasing functions of the asymptotic 
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purities and decreasing functions of the couplings to the baths. This should 
be expected, recalling the well understood synergy between entanglement and 
purity for general quantum states: the ideal vacuum environment, whose deco- 
herent action is entirely due to losses, is the one which better preserve both the 
global information of a state and its correlations. 

As we have seen, two-mode squeezed thermal states constitute a relevant 
class of Gaussian states, parametrized by their purity fi and by the squeezing 
parameter r according to Eqs. (|59|l . In particular, two-mode squeezed vacuum 
states (or twin-beams), which can be defined as squeezed thermal states with 
/i = 1, correspond to maximally entangled symmetric states for fixed marginal 
purities |17j . Therefore, they constitute a crucial resource for quantum informa- 
tion processing in the continuous variable scenario. For squeezed thermal states 
(chosen as initial conditions in the channel), it can be shown analytically that 
the partially transposed symplectic eigenvalue i/_ is at any time an increasing 
function of the bath squeezing angle ip2'- "parallel" squeezing in the two chan- 
nels optimizes the preservation of entanglement. Both in the instance of two 
equal squeezed baths {i.e. with ri = r2 = r) and of a thermal bath joined to a 
squeezed one (i.e. ri = r and r2 = 0), it can be shown that is an increasing 
function of r ISU]. Such analytical considerations, supported by a broader nu- 
merical analysis, clearly show that a local squeezing of the environment faster 
degrades the entanglement of the initial state. The same behavior occurs for 
purity. 

In order to illustrate the behaviour of the nonclassical depth r in the noisy 
channel, let us consider standard form states evolving in thermal environments. 
For simplicity, let us assume 71 = 72 = 7. According to Eqs. (|16() and 1)6 5() . one 
has, for the evolving nonclassicality (recalling that |c2| > |ci|) 



r{t) =l-ka + b) c^^' - /^2fL±i^(i _ e"-*) 




(a - 6) e-7* + (1 - e-^*)) + 4c2 e-27t 



(70) 



This function is a decreasing function of the parameters ^i^o ■ the thermal noise 
contributes to destroy the nonclassical features of the initial state. To study 
the effect of the squeezing of the bath on the nonclassical depth, we specialize 
to the instance of two mode squeezed thermal states, which are an archetypical 
class of nonclassical two mode states, characterized by squeezing in combined 
quadratures. In this case it can be easily shown that, in order to minimize the 
smaller eigenvalue of cr (thus maximizing r) the choice (^2 = is optimal. We 
will thus make such a choice in the following. The nonclassical depth of the 
initial two mode squeezed state ^^i^^ in a channel with parameters fii. 
for i = 1,2, takes the following form 

^1 _ cosh(2r) _ e'^-^-^i2^+ e-^--'-^^^^ _ 

2 2y/JI Vloo/i2oo 




= /i2oo-e ''•^-Mioo ' ^ sinh(2r)2 



2/iloo/^2oo 



(71) 
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Figure 13: Evolution of the nonclassical depth of an initial two-mode squeezed thermal state 
with fi = 0.9 and r = 1.5. The dashed line refers to the evolution in a non squeezed bath with 
MicxD = M2oo = 0.5; the dot-dashed line refers to a non squeezed bath with /iioo = 0.25 and 
M2oo = 1; the dotted line refers to a bath with /^loo = M2oo = 0.5 and ri = r2 = 0.2; finally, the 
continuous line refers to a bath with fj-iao = M2oo = 0.5, rioo=o.6 3nd r2oo = 0. 



Eq. H71(l reduces to the following simple form for the evolution in equal baths 
(with ^loo = /i2oo = ^/^^ and rioo = ^200 = ?'oo) 

1 - c-^* - 21!:^ fi _ e-rt) 



As can be seen in Fig. ^1 the local squeezing of the baths, reducing the quan- 
tum noise in one quadrature of the multimode system, drastically increases the 
duration of the nonclassicality of the state and, generally, the value of its non- 
classical depth at any given time. This is due to the symmetry of two mode 
squeezed states under mode exchange: such states can take advantage of reduced 
fluctuations of any quadrature of the bath. Interestingly, while the nonclassical 
depth is enhanced by the local squeezing of a quadrature (thus implying an 
improved preservation of nonclassical features like subpoissonian photon num- 
ber distributions), the entanglement is not. This is due to the intrinsically non 
local nature of the entanglement: the advantage which could be achieved by 
squeezing a local quadrature is balanced by the increased fluctuations in the 
conjugated quadrature, which usually makes squeezing not favourable to the 
aim of preserving entanglement. 



8 Concluding remarks 

We have carried out a quantitative analysis of decohercncc of continuous variable 
systems interacting with general Gaussian environments and reviewed many re- 
lated results. The method we have presented to study the decohercncc rate may 
be applied to other systems of interest, like qubit systems under non unitary 
evolutions. Several relevant configurations have been considered and exhaus- 
tively analysed, characterizing their rate of decoherence by keeping track of the 
decay of the global degree of purity, of indicators of nonclassicality and, for two 
mode states, of quantum and total correlations. 
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Quite in general, we have shown that, as long as one restricts to the Gaussian 
setting, squeezing the bath (or, equivalently, the intial state while letting the 
bath being thermal) does not help to better preserve either the overall coherence 
of the state or its quantum correlations. However, such a squeezing proves 
effective in delaying the decoherence of more deeply non classical states, like 
cat-like states resulting from coherent superpositions of Gaussian states or of 
number states. Furthermore, quite interestingly, wc have shown that a local 
squeezing of the baths may improve the preservation of the mutual information 
in two-mode systems. 

We remark that our results are of direct interest to recent developments in 
experimental quantum optics, especially related to quantum information and 
quantum control. Indeed, a crucial step towards the development of quantum 
information technology is the achievement of a sufficent quantum control capa- 
bility, i.e. of the ability of engineering quantum signals and feedback techniques 
acting on the dynamics of a quantum system. In fact, the implementation of any 
quantum information protocol relies on maintaining quantum coherence in the 
system for a significant period of time and so requires some kind of mechanism 
to eliminate or mitigate the undesirable effects of decoherence. In this frame- 
work, a precise knowledge of the decoherence dynamics is desirable, especially 
in the continuous variable regime, where the field of quantum control originated 
and has a strong experimental impact |82l l8Hj . 

In order to make this point clearer, let us explicitly regard the following ex- 
ample. Consider the continuous variable teleportation of a single-mode coherent 
state by exploiting a two-mode squeezed thermal state as entangled resource 
(for a detailed description of the protocol, see Ref. |H11)- Now, it may be shown 
|SH] that the optimal teleportation fidelity F (averaged over the whole complex 
plane) for such a protocol is given by a simple function of the smallest partially 
transposed symplectic eigenvalue of the two-mode squeezed state: 

If the two modes which share the entangled state are, say, stored in two distant 
cavities, waiting to be used, the decoherence they experience will gradually cor- 
rupt the fidelity of the teleportation protocol. Our study allows to keep track of 
the quantity i>- during the dissipative evolution of the state as a function of var- 
ious environmental parameters, and thus to exactly determine the teleportation 
fidelity achievable as a function of time. For instance, considering an initially 
pure shared two-mode squeezed vacuum with squeezing parameter r, evolving 
in two environments with, for simplicity, the same coupling 7 and asymptotic 
purity ^00, one gets 

= 1 + C-^^r-.t + (1 _ e-7t)/^^ ■ 

Notice that such a result takes into account both losses and thermal noise. In 
the more general instance, let us remark that the entanglement time, extensively 
analyzed in Sec. [see Eqs. (|67I69|) ] and which may be analytically determined 
following the approach we have presented, coincides with the time over which 
quantum teleportation allows to beat the classical fidelity, equal to 0.5, as shown 
by Eq. (|73|l (at t^nt reaches 1/2 and then keeps increasing). After such a 
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time the entanglement is gone because of local decoherence: the shared resource 
becomes useless to quantum informational aims. 



A Determination of mixedness and entanglement 
of two-mode states 

Here we provide explicit expressions which allow to determine the exact evo- 
lution in uncorrelated channels with 71 = 72 = 7 of a generic initial state in 
standard form. The relevant quantities £V, 1^, Sv, I and r are all functions 
of the four S'p(2,r) © Sp(2,R) invariants Dcta. Dct/3 Det7 and Deter. Let us 
then write these quantities as follows 



Det o- = Efe c"''^^* , (75) 

fe=0 

2 

Beta = ^afcC"''", (76) 



fe=0 

2 



Det/3 = ^/Jfec-'^"^ (77) 
Det 7 = 72 6-2", (7g) 
defining the sets of coefficients E^, a^, Pi, 7^. One has 

2,2 2jCOsh2r2 2COsh2ri cosh 2ri cosh 2r2 cosh2ri cosh2r2 

a b + + —J — a b ab h ab a— 5 

^2 M2 Ml M1M2 4/ilA*2 VlAi2 

2 9x / cosh2r2 &cosh2ri cosh 2ri cosh 2r2 sinh 2ri sinh 2r2 cos 2o32 
+ + C2) a — \ ab 

V 2^2 2/ii Vi/i2 Vl/i2 

n .J / sinh 2r2 cos 2c52 ,sinh2ri sinh2ri cosh2r2 cosh 2ri sinh2r2 cos 2g92 
+(c-, — c,) a h — 

V 2/X2 2/^1 4/ii^2 ViAi2 



" 4Ai2 4/^2 /i2 



- 2- 



2 cosh 2ri cosh 2ri cosh 2r2 „ cosh 2ri 
ao 2ab \~ 3a- 



36 



Ml 



M1M2 



cosh 2r2 
4MfAi2 



sinh 2r2 cos 2(^2 sinh 2ri sinh 2ri cosh 2r2 
a h — 2 ^ 2 



2/i2 



2^1 



4miM2 



cosh 2ri sinh 2r2 cos 2iy92 
4miM2 



-(c! +c^) ( a 
1 



cosh 2r2 6 cosh 2ri cosh 2ri cosh 2r2 sinh 2ri sinh 2r2 cos 2iy92 



2^2 



2fii 



4miM2 



4miM2 



(80) 
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cosh 2ri cosh 2r2 cosh 2ri cosh 2r2 
+ "T^ + ao Sa— 5 36— -^5 

, n n, / cosh 2ri cosh 2r2 sinh 2ri sinh 2r2 cos 2o92 



4/iiAi2 4/ii^2 
n n / sinh 2ri cosh 2r2 cosh 2ri sinh 2r2 cos 2o92 \ 1 

' 4^i^2 4/ii/i2 / 16^f^^ 



cosh2ri ^cosh2r2 
1 

16m?M2 



So — 2 2 ' ^^^) 



2 cosh2ri 1 , ^, 

a2 = a-^ - a h — , (84) 



Ml 4^f 
cosh2ri 1 
Ml 4/x2 

1 



ai = a 1-2—, (85) 



"0 = 77^ , (86) 



, 9 , cosh 2r2 1 
/32 = 6' - 6 ^ + -3 , (87) 



A'2 4^2 
cosh2r2 1 
M2 4/x, 



A = .. -2,2 



/3o = ^, (89) 

72 = C1C2 . (90) 

The coefficients of Eq. (jfi6|l . whose solution fcent allows to determine the 
entanglement time of an arbitrary two-mode Gaussian state, read 

(91) 
(92) 

I72I , (93) 

(94) 

i • (95) 
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